Photon Statistics for Single Molecule Non-Linear Spectroscopy 
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We develop the theory of non linear spectroscopy for a single molecule undergoing stochastic 
dynamics and interacting with a sequence of two laser pulses. We find general expressions for the 
photon counting statistics and the exact solution to the problem for the Kubo- Anderson process. In 
the limit of impulsive pulses the information on the photon statistics is contained in the molecule's 
dipole correlation function. The selective limit, the semi-classical approximation, and the fast 
modulation limit, exhibit rich general behaviors of this new type of spectroscopy. We show how the 
design of external fields leads to insights on ultra-fast dynamics of individual molecule's which are 
different than those found for an ensemble. 

PACS numbers: 82.37-j, 82.53-k, 05.10.Gg, 42.50.Ar 



Nonlinear optical interaction of a sequence of laser 
pulses with matter provides a powerful tool for the in- 
vestigation of dynamics of ensembles of molecules in a 
wide variety of chemical, physical and biological sys- 
tems Recently, van Dijk et al [3 reported of the 
first experimental study of an ultra-fast pump-probe sin- 
gle molecule system. Unlike the previous approaches to 
non-linear spectroscopy where only the ensemble average 
response to the external fields is resolved, the new ap- 
proach yields direct information on single molecule dy- 
namics, gained through the analysis of photon counting 
statistics. Although the original experiment 0] was con- 
ducted on a molecule undergoing a relatively simple re- 
laxation process, the huge potential of combining non- 
linear spectroscopy with single molecule spectroscopy, in- 
spires many unanswered theoretical questions: What are 
the fundamental physical limitations of the investigation 
of fast dynamics when spectral selectivity (defined be- 
low) is limited? How does the information contained in 
these experiments differ from the information contained 
in simpler continuous wave experiments? What is the 
finger print of coherence in these type of experiments, 
and precisely how its infiuence on photon statistics is 
suppressed due to dephasing processes? Finally, how to 
design the external control fields so that new information 
on dynamics of molecules is gained. While the answer 
to these questions, depends on the particular dynamics 
of the molecule under investigation, we present a theory 
based on the Kubo- Anderson model, which yields general 
insights on the problem. 

Consider a sequence of two laser pulses interacting with 
a single quantum system as a molecule, an atom, or a 
nano-crystal. The pulses are assumed to be very short 
compared with the radiative life-time of the emitter, so 
that the probability of photon emission during a pulse is 
negligible, hence a pair of pulses yields two photons at 
most. Repeating the experiment many times one may 
obtain the probabilities Po,Pi and P2 of emitting 0,1 
and 2 photons. In this manuscript we investigate the re- 
lation between Pq , Pi and P2 and the dynamics of the 
underlying system interacting with the external fields. 



In turn this type of photon statistics reveals important 
information on fast dynamics of molecules in the con- 
densed phase, information which is very difficult to ob- 
tain using other theoretical approaches to single molecule 
spectroscopy [1, 0, S H 0) Q ■ We consider a model of 
a single molecule undergoing spectral diffusion and in- 
teracting with a pump-probe set up and show that de- 
pending on characteristics of the stochastic process and 
laser field parameters, different types of nonlinear spec- 
troscopies emerge. In particular, sensitivity to the phase 
accumulated by the molecule in the delay period is found 
and impulsive and selective type of spectroscopies are 
considered in detail. Spectral diffusion used here as a 
prototype of the molecule's dynamics, is found in many 
molecular systems 0, 0, [13, El , is easy to detect using 
the spectral trail technique, when the process is slow. 
Our goal is developing general methods suitable for the 
detection of a wider range of dynamics. 

The treatment begins with a two level molecule inter- 
acting with a time dependent laser field according to the 
optical Bloch equation [l[ 



cr = L{t)a + t(7 



(1) 



where a — (cco, Cgg, CTgc, Cog)- As usual Coo and Cgg are 
the populations of the excited and ground states, and 
f go , CTog are the off diagonal matrix elements of the den- 
sity matrix. The operator 
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describes the interaction of the molecule with the driving 
electromagnetic field through ri/(i) where f2 is the Rabi 
frequency, while io(€) is the time dependent absorption 
frequency of the molecule soon to be discussed. Finally 
the operator F = r|g)(e| describes the transition from the 
excited state into the ground state, due to spontaneous 
emission with F designating the emission rate. 

The spectral diffusion process ijj(t) is modeled using the 
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Kubo- Anderson approach: uj{t) = u}o + 5w{t) where ujq is 
the bare absorption frequencyof the molecule, and Sw(t) 
is a random function of time [H, • We will assume that 
this process is stationary, its mean is zero, its correlation 
function is {Sw{ta + t)Sw(to)) = v'^%j}{t) with V'(O) = 1 
and -0(00) = 0. We will later demonstrate our results 
using the Kubo- Anderson process 0, where uj{t) — 
uiQ + v or Lo{t) = loq — V, with the rate R determining the 
transitions between the + and — states, which was used 
to model single molecules in low temperature glasses [^. 

The photon statistics is obtained from the solution of 
Eq. ([1]) in terms of an iterative expansion in F Q . This 
yields a{t) — X!J^o ^^"H^) where 
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The Green function g{t, t') = f exp L{t)dt [f - the 

time ordering operator) is the evolution operator in the 
absence of the spontaneous emission (i.e. without F) and 
therefore cr*^"^ {t) describes the state of the system condi- 
tioned by n emission events. The natural basis we shall 
use is |e) = (1, 0, 0, 0) which means that the molecule is in 
the electronic excited state, |g) = (0,1,0,0) the ground 
state, while |c) = (0,0,1,0) and |c*) = (0,0,0,1) de- 



scribe the coherences. Hence, the probability of emission 
of n photons is P„(t) = (elCT^"^*)) + (gk^"^^), which 
we calculate below in the limit of long measurement time 
t 00. 

In our model two identical pulses: pump and probe 
interact with the molecule, 
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where ujl is the laser frequency, ti = — t2 — T are 
pulses duration, and A = i2 ^ is the delay between 
the pulses. Since the Green function can be written 
as a product of propagators of shorter time intervals: 

g{t,o) = g{t,h)g{h,t2)g{t2,ti)g{h,o), it is possible to 

treat everyone of them separately by dividing the range 
of integration over time in Eq. ([31), and inserting closure 
relations in appropriate places. When the laser is turned 
off the Green functions g{t2,ti) and g{t,t3) are nearly 
trivial, and one may find general expressions for P„ (any 
n) for a particular realization of the spectral diffusion 
process. In the limit of two short pulses obeying FT <C 1 
we neglect the processes where photons are emitted dur- 
ing the pulse events, and find 



P„ = P^'^ [CO {h),u;{t2)] + {e''^A-rA/2^Coh (i^) (i,)] + c.C.} , (5) 

here 9 A — J^^^'^'^ uj{t)dt is the random phase accumulated in the delay interval. Pn^^ given in Table 1 be- 
low, describe the semi-classical paths in the sense that they include only the transitions between the pure ex- 
cited and ground states, whereas represent the contribution of the coherence to the photon statistics. 
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Table 1 

r 



In our further calculations we assume that changes in 
the absorption frequency of the molecule during the pulse 
events are negligible, namely that the rate of the spectral 
diffusion process changes R satisfies RT <C 1. As a conse- 
quence the random absorption frequency of the molecule 
at the moment of excitation will be taken as uj(t) = Lo{ti) 
during the first pulse event and uit) = Lu{t2) during the 



second pulse. If this condition is met, the Green functions 
needed for the calculation of the matrix elements in Ta- 
ble 1 are found within the rotating wave approximation 

(RWA) gHh)] = g{h,0) - e^p{L^'^^[u;{h)]T} 
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and similarly for Q [^(^2)], 
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where S{ti) = lvl — is the detuning. Otherwise, the 
pulses yield time average information on the dynamics. 

The results in Table 1 describe the possible physical 
paths of photon emission. For example consider the first 
term of P^^^ in Table 1, the molecule starts in the elec- 
tronic ground state |g), it then evolves with the Green 
fimction of the first pulse Q [ujlti)] without emitting a 
photon to the excited state |e), it remains in the excited 
state in the delay period without emitting a photon [with 
probability exp(— FA)], and then the second pulse stim- 
ulates the molecule from |e) to |g). From Eq. ([5]) we 
see that the non-classical terms are important only 

when FA is not too large, as expected. Also Aj"'' = 
in Table 1, since the emission of two photons from two 
short pulses is possible only when one photon is sponta- 
neously emitted in the time interval between the pulses 
and the second photon after the second pulse, hence in 
this case the coherence is lost due to the collapse of the 
wave function. 

When FT — > normalization condition is Pq + 
Pi + P2 = 1 as expected. In this case the sponta- 
neous emission does not contribute during pulse events, 
and reversibility and symmetry of matrix elements is 
found: {e\g[Lo{h)]\c) = (g|^[^(ti)]|g), {e\gHh)]\g) - 



(g|5[(jj(<i)]|e) etc. It is then easy to use Table 1 
and show that Po^'^[w(ti), ^(^2)] + PP''[w(ti), ^(^2)] + 
P2^^[uj{ti),uj{t2)] = 1, SO the classical paths conserve 
probability, and Ai°^ = —Aq°^. These simple relations 
are valid for any realization of the underlying stochastic 
process and they are valuable in reducing the number of 
matrix elements needed for the calculation of the photon 
statistics for a specific stochastic path. 

To obtain the solution of the problem one must take av- 
erages of Pn in Eq. ([5]) over the stochastic process. This 
requires knowledge of the joint probability density func- 
tion (PDF) P[u}{ti),u;{t2),9A] of finding the molecule's 
absorption frequency in the infinitesimal range near uj(ti) 
at ti, near ^(^2) at ^2 with accumulated random phase 
OA. Before considering the concrete example of the Kubo 
Anderson process, we consider different limits allowing 
us to make a number of conclusions valid for all types of 
spectral diffusion processes. 

Impulsive limit <C J7. - Taking the limit T — > 0, 
f2 cxD in the pulse Green function Eq. ([6]) while f2T 
remains finite in such a way that v <ti D, we define the 
impulsive limit. The matrix elements of and 
g[ijj{t2)] become independent of the value of Lu{t) at the 
moment of the excitation, provided that the laser de- 
tuning is small compared with Q. Thus instead of the 
multi variable PDF P[LLj{ti),uj{t2),0A] we now have to 
deal only with the one variable PDF of the phase 9 A. 
As a result the photon statistics shows an interesting re- 
lation with linear continuous wave spectroscopy. Using 
Eqs. ([5l[6]) and Table 1 we find for stationary processes 
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where 

0(A) = (exp(i / buj{t!)At!) (8) 
Jo 

is the well investigated Kubo- Anderson correlation func- 
tion whose Fourier transform is the line shape of the 
molecule according to the Wiener-Khintchine theorem 
[1^ . This type of impulsive limit offers no spectral res- 
olution in the sense that the photon statistics Eq. ^ is 
sensitive only to the random phase 9 accumulated in the 
time interval between the pulses and not to the temporal 
state of the molecule at the time of the pulses. From Eq. 



([7]) we see that for a 7r/2 pulse with SIT — tt/2 the impor- 
tance of the coherence terms and hence the correlation 
function 0(A), on the photon statistics is the strongest, 
since the n/2 pulse excites the off diagonal terms. 

Semi- classical approximation. - The influence of coher- 
ence on photon statistics in many experimental cases is 
expected to be difficult to detect. Either due to de- 
phasing effects or simply because of the large value of 
the optical transition frequency. In these cases a practi- 
cal approximation is to remove all the coherence terms, 
and leave only the semi-classical terms P^'**. From Ta- 
ble 1 it is easy to see that in this case we must deal 
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with the two dimensional PDF P[u!{ti), ^(^2)], instead of 
P[uj{ti) , U!{t2) , 0^] for the general case. And then 

/•OO /"OO 

(P„) = (/^f ^> = / / Pf^PKii),c.(i2)]da;(ti)d^(t2). 
Jo Jo 

(9) 

Below we give an explicit example for this case for the 
two state Kubo- Anderson model. 

Exact solution. - Now we obtain the exact solution for 
the two state Kubo- Anderson Poissonian process, where 
the absorption frequency of the molecule jumps between 
a + and — states where uj = loq ± v. We denote the 



initial state, during the first pulse with i = + or j = — , 
similarly the final state of the molecule at the second 
pulse is / = + or — . Since the random phase A6' = wqA-I- 
i^(r+ — r~) where are occupation times in states + 
and - [13], satisfying A = T+ + T' , the joint PDF 
P{iJ,A9) is equivalent to the joint PDF h{iJ,T+) of 
finding the molecule in state i in the first pulse, state / in 
the second and with the occupation time T+ between the 
two pulses. We leave technical details on the calculation 
of h{i,f,T~^) for a longer publication, taking averages 
over Eq. ^ gives: 



«=±,/=± 

where h {i, /, ~2ih') is the Laplace T+ 2iiy transform of h{i, /, 
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while V^^ = [1 + exp(-2i?A)]/4, = [1 - 

exp(— 2i?A)]/4, are probabilities of finding the particle 
initially in state i and finally in state /, which are easy 
to obtain from Poissonian statistics. 

As an illustration we apply Eq. (|10p to the case of two 
infinitely strong UT — n pulses, so no coherence is built. 
At first assume that f2 <C and ujl = ujq + 1' so that the 
laser is in resonance with the + state and not with the 
— state. With this selective limit unlike the impulsive 
limit, temporal resolution is found 
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(12) 

and (P2) = 1 — (Pq) — {Pi)- These results make perfect 
physical sense, because the molecule emits a single pho- 
ton only when it is once in state -I- and once in state — , 
hence (Pi) = +V '". These simple solutions show 
that for weak fields ^ i/ the photon statistics does 
not depend on frequency shifts v, and hence to explore 
the dynamics of the molecule we must consider stronger 
fields. Using Eq. pO)) for tt pulses, we find the exact 
expression 
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which reduces to Eq. in the limit ^ <^ v. In the 

opposite limit ^ v we have (Pi) = 0. Hence for the 
investigation of spectral shifts and rates we cannot use 
neither too weak or too strong fields which with the con- 
dition QT = TT means that pulses must not be chosen 
arbitrarily short. Eq. (fT5|) shows precisely what fields 
yield information on the process beyond simple limits. 

Fast modulation. - An interesting case is the fast mod- 
ulation limit R ^ v, where motional narrowing effects 
take place. We consider the limit i? — > 00 and ^ 00, 
in such a way the the spectral diffusion dephasing rate 
FsD — jR is finite. The latter is a measurable physi- 
cal observable since it gives the width of the line shape in 
continuous wave spectroscopy [l^ . To obtain the photon 
statistics in this limit we can proceed by expanding the 
exact solution Eq. (|10p in terms of the large parameters 
R and v. However, there is a more general and sim- 
pler approach. Notice that v <^ R <^ 1/T and to have 
a finite probability of photon emission must remain 
a constant of order unity, therefore u <^ Hence in 
this fast modulation limit the pulses must be impulsive. 
Then Eq. Q holds with 4>{A) = exp(-rsD A/2), which 
means that (Fgo + r)/2 is the renormalized decay rate 
which damps the coherence terms. This result is valid 
for many types of fast spectral diffusion processes, and is 
not limited to the exactly solvable two state process. 

Semi- classical selective limit.- Neglecting the coher- 
ence terms and taking a laser in resonance with the + 
state and out of resonance with the — state i.e. ojl = 
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+ ^ we find from the exact result Eq. (fTO|) Acknowledgment This work was supported by the 

Israel Science Foundation. 

(Po) = cos2 j +sin4 j (1 + e-^^) ^ , 



(Pi)=sm^(-j-sm^(^-j^^ (14) 

and (P2) = 1 - (Pi) - (Po)- Note that these results 
exhibit Rabi oscillations and that (Pi) is independent of 
exp(— FA) the latter behavior is general: for any spectral 
diffusing process (Pf^'*^) is independent of exp(— FA) in 
the limit of short and strong pulses due to the mentioned 
symmetries of the matrix elements of the pulse Green 
function. 

Summarizing perspectives.- Theoretical investigation 
of the new field of single molecule non- linear spectroscopy 
was presented. We showed that in the limit of impulsive 
pulses information on the photon statistics is given by 
the Kubo- Anderson correlation function, for any station- 
ary spectral diffusion process. To obtain information on 
spectral shifts and spectral rates beyond the trivial lim- 
its, one must choose carefully the duration of pulse and 
its strength, as our analytical solutions clearly demon- 
strated. Our results can be checked in experiments in 
low temperature glasses, and our methods provide the 
theoretical basis for the investigation of fast dynamics 
on the single molecule level, thus opening the door for a 
vast number of applications. 
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